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ABSTRACT 

£> | We present the mid- infrared spectra of seven of the most powerful radio-galaxies known 

C*~) . to exist at 1.5 < z < 2.6. The radio emission of these sources is dominated by the 

AGN with 500 MHz luminosities in the range 10 27 8 -10 29 1 WHz- 1 . The AGN signa- 
ture is clearly evident in the mid-infrared spectra, however, we also detect polycyclic 
aromatic hydrocarbons emission, indicative of prodigious star formation at a rate of up 
to ~ 1000 M yr . Interestingly, we observe no significant correlation between AGN 
power and star formation in the host galaxy. We also find most of these radio galaxies 
CNI ■ to have weak 9.7 /im silicate absorption features (rg.7 Mm < 0.8) which implies that 

their mid-IR obscuration is predominantly due to the dusty torus that surrounds the 
central engine, rather than the host galaxy. The tori are likely to have an inhomo- 
geneous distribution with the obscuring structure consisting of individual clouds. We 
estimate that these radio galaxies have already formed the bulk of their stellar mass 
and appear to lie at a stage in their evolution where the obscured AGN dominates the 
energy output of the system but star formation is also prevalent. 

Key words: galaxies: high-redshift, active, nuclei, starburst, quasars: general 
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1 INTRODUCTION 

It is now understood that supermassive black holes 
(SM BHs) are at the c e ntres of all massive galax- 
ies |Magorrian et all Il998l : iRichstone et all Il998j). These 
SMBHs can be activel y accreting (B landfor d &: Pavndfl982l ; 
IChurazov et"afl 120051 : iRussell fc Fender! l2010h which has 
given rise to the populatio n of active galacti c nuclei (AGN). 
A small fraction (< 10%, iTasse et al.ll2008l ) of these AGN 
are 'radio-loud' in that their radio flux is greater than their 
i?-band flux by a factor > 10. More specifically however, 

* E-mail: jir28mssl .ucl . ac .uk 



the radio- loud frac tion may vary wi th redshift, mass and 
optical luminosity (I Jiang et al.ll2007l ). Powerful radio- loud 
AGN show the presence of twin jets of relativistic parti- 
cles which emit syn chrotron radiation at radio wavelengths 
l|Lvnden-Belllll969l '). The radio jets can reach up to hun- 
dreds of kilo-parsecs in length and extend well beyond the 
realms of the host ga l axy. Under the unification model 
( Rowan- Robinson! [l977l : lLawrenceHl987l h radio galaxies are 
the hosts of radio-loud type-2 AGN, whose jets are oriented 
approximately perpendicular to our line of sight a nd much of 
the central engine i s obscured by a dusty torus (|Antonuccil 
1984; Barthcl 1989), allowing the host galaxy to outshine 
the AGN in the optical. 
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In the local Universe, radio galaxies are associated 
with massive giant elliptical galaxies (Mat thews et al.l 

1 1964). Studies of r a dio galaxie s at high redshift 

jPentericci et all l200l] ; IZirm et al.l 120031 ; ISevmour et~ai1 



120071 ; |Pe Breuck et al 



indicate that they too are 
massive systems (10 10 -10 12 Mq) up to at least z ~ 5.2. 
As a result, it has been suggested that high-redshift radio 
galaxies (HzRGs) are the precursors of the modern day 



llArchibald et aU 


200 ll; Stevens et al. 20031; Reuland et all 


2004; iMilev & De Breuck 20081; Galametz et al. 20121; 


Mavo et al.ll2012f) 


. HzRGs have been found to host SMBHs 


of 10 9 -10 10 M Q 
higher accretion 


|Nesvadba et al.l l201ll). Thev also have 
rates than their local counterparts as 



inferred from mid-infrared (mid-IR) continuum luminosities 



e.g. rest-frame vL v (15 /tm) > 10 erg s 



Ogle et al. 20061; see als o Seymour et al 



2007al ; lDe Breuck et aLlkoid ; I Ogle et"^Tll200^ . T2OI0I T 



f or z >0.5, 
Saiina et al.l 



From other investigations into the host properties of 
HzRGs, there is tentative evidence that some have very 
high star formation rates (SFRs). For example, from sub- 



mm observations of r adio galaxies jDunlop et al. 19941: 



Archibald et all l200ll ; iReuland et alj |2004 iGreve et al.l 



have been reported, while 
derived from the lu- 



20061 ). SFRs of ~ 1000 M yr" 
SFRs of up to 1500M Q yr _1 were 
mi nous (>100kpc) Lyq halos around some ra dio galax- 
ies llReuland et al.ll2003j;IVillar-Martfn et alj|2003l ). Further- 
more. ISevmour et al.l (|201lT ) investigated the star-forming 
properties of radio-loud AGN in the far-IR and while they 
made no distinction between type-1 and type-2 AGN, for 
their high-redshift (1.2 < z < 3) sample they found a 
mean SFR range of 80-580 Mq yr _1 . 

The spectral energy distributions (SEDs) of radio 
galaxies show copious amounts of emission from both the 
AGN and host galaxy. In the mid-IR, along with re- 
process ed emission from dust associated with the AGN 
torus l|Neugebauer et all Il987l) . pol ycyclic aromatic hy- 
drocarbon (PAH) emissi on features (jPuget fc Leger!ll989l ; 
lAllamandola et all D-999) may be present. PAH emission 
manifests itself as broad (~ 1 pm) features, primarily at 6.2, 
7.7, 8.6 and 11.3 /im, that pr edominantly appear in the mid- 
IR sp ectra of dusty galaxies (|Helou et al.ll2000l ; ISmith et al.l 
2007). The various studies on AGN-starburst composite 
sources show that PA Hs can survive in galaxies that con- 
tain an active nucl eus (|Desai et alj|2007l ; lArmus et al.ll2007i ; 
iTrever et al.l I2OI0I ). PAH emission can be used as an ap- 
proximate measure of a galaxy's SFR as it is thought 
that PAHs become excited by the absorption of UV emis- 
sion from massive, young, hot stars. The relation is fur- 
ther supported by the observed correlation between PAHs 
and the cold dust reservoirs associated with star forma- 
tion (lHaas et al.ll2q02l; iPeeters et alj|2004 lLacv et al.ll2007l; 
ISaiina et al.l |200S| ; iPope et allbooa hereafter P08). It has 
been sug gested that PAHs could b e present in local radio 
galaxies (jSiebenmorgen et al.l |2004| ) . Ho wever, their weak 
detection or ab s ence i n other studies dHaas et alj 120051 ; 
lOgle et all 120071 . |2010| ; iNesvadba et~aH |2010| ) suggest con- 
siderably lower SFRs (<5M Q yr _1 ) than those detected 
in radio-loud AGN at hi gh-redshift (ISaiina et al.l l2007bl ; 
ISevmour et~ai1l200Sl . l201ll ; lOgle et al.ll2012D . From the few 
mid-IR spectra of HzRGs to date, P AH emission ha s been 
detected on a couple of occasions l|Sevmour et al.l 120081 ; 



lOgle et all |2012| ). This has given SFRs of up to a few 
1000 M©yr -1 which compliment the sub-mm based esti- 
mates. 

Another feature that may be present in the mid- 
IR spectra of radio galaxies is silica t e ab s orption at 
9.7 /nm and 18 /mi dGenzel fc Cesarsksi |2000| ; lYan et al 
20071 ; ISaiina et ail l2007al ; iNikutta et alj |2009| ; ISales et al 



20111 ). From the prevalence of weak silicate absorption 
features observed in type-2 AGN, it has been suggested 
that their obscuring tori are characterise d by a clumpy, 



rather than smooth, di s tribution of dust dLevenson et al 



2007 ; ISpoon et alj [2007 ; Im anishi et al. 2007; ISaiina et al 



2007a). Strong absorption features on the other hand, 



have been interpreted as evidence for a large contribu- 
tion f rom the host ga laxy to the obscuration of the nu- 
cleus (|Ogle et al.ll2006l ). The weak silicate absorption fea- 
ture at 9.7 /im for one HzRG; 4C 23.56 l|Sevmour et all 
2008], T9.7 0.3 ± 0.05) is comparable to those found 
for some lower-redshift and less- luminous radio galaxies 
l|Qgle et al . 2006; Clcary et al. 2007). However, a small num- 
ber of sources from lOgle et alj l|2006l ) showed deeper silicate 
absorption as did t he high-redshift radio-loud AGN from 
ISaiina et al.l (|2007bi . rg.7 Mm > 1). However, the latter sam- 
ple was selected to contain sources with significant obscura- 
tion. 

In this paper, we present the mid-IR spectra 
(5/xm < A rcst < 13pim) of seven HzRGs, the largest sam- 
ple to date, t hat were obtained using the Infrared Spec- 
trograph (IRS. lHouck et al.|[ 2004) instrument on-board the 
Spitzer Space Telescove (|Werner et al.l |2004 ). We investi- 
gate the connection between star formation and the radio- 
loud phase of powerful obscured AGN by decomposing the 
spectra into AGN and host galaxy components. This al- 
lows the relative AGN and star-formi ng contributions to 
the mid-IR output to be determ i ned ijSaiina et alj l2007al : 
P08; [ Menendez-Delmestre et al 
Coppin et alj 2010; Petric et al 



20091, hereafter MD09 
I2OIII ; iFiolet et al.1 l20ld : 



Mullanev et alj|201ll ). 

The paper is structured as follows: In Sj2]we describe the 
sample, the observations and our data reduction pipeline, in 
33 we present our spectral fitting model, while in Uand 33 
we comment on our results and draw conclusions. 

Throughout, we assume a A-dominated flat Universe, 
with H = 70kms" 1 Mpc" 1 , « A = 0.7 and fi m = 0.3. Unless 
otherwise stated, uncertainties are 1 a and upper limits are 
3a. 



2 THE SAMPLE, OBSERVATIONS AND 
REDUCTION 

2.1 Sample selection 

Our sample of radio galaxies originate from the 
Spitzer High-r e dshift Radio Galaxy ( SHzRG ) project 
|Sevmour et all 120071 ; |Pe Breuck et al.1 l2010l . hereafter 
DB10) which focused on targeting the host galaxy prop- 
erties of some of the most powerful obscured radio-loud 
AGN in the Universe. The SHzRG sources were selected 
fro m multiple radio gala xy surveys (for list and references, 
see lSevmour et al"1l2007l ) for which radio flux densities were 
obtained from all-sky, low-frequency (74, 327/352, 365 and 
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1400 MHz) surveys jRengelink et aL 1997; Douglas ct al. 
1 19961 ; ICondon et all 1 19981 ; ICohen et all l2007f ). They were 
selected to uniformly sample the redshift-radio luminosity 
plane for z > 1 and L3 ghz > 10 26 WHz -1 . Radio galaxies 
with ancillary data were preferred. With these selection cri- 
teria, this gave a sample of 70 objects that is representative 
of the overall HzRG population. Spectroscopic redshifts and 
rest-frame 3 GHz luminosities were in the range 1 < z < 5.2 
and L 3 GHz = 10 26 ' 2 -10 28 ' 7 WHz" 1 , respectively. Due to 
their extremely high radio luminosities, objects of this 
nature are very rare, with only a few hundred known at 
2 > 1 acr oss the whole s ky. Most have steep spectra in 
the radio |Ker et al.l |2012| ) with th eir 3 GHz synchrotro n 
emission free of Doppler-boosting l|Sevmour et alj l2007h . 
Rest-frame 500 MHz luminosities, Lsoomhz, were calculated 
by interpolating between the luminosities from the differ- 
ent surveys. The isotropy of the synchrotr on emission is 
expe cted to anti-correlate with frequency l|Blundell et al.l 
1998) and 500 MHz was the lowest rest-frequency, for these 
redshifts, at which luminosities could be interpolated. The 
SHzRG objects were first observed with all seven photo- 
metric bands of t he Spifzer-Infrared Array Camera (IRAC, 
iFazio et all |2004|)_ and Multiband Imaging Photometer for 
Spitzer (MIPS. iRieke et al]|2004h instruments. Later, for a 
select few, we used the IRS for follow-up spectroscopy. 



The sample presented in this paper consists of seven 
of the eight SHzRG objects that have been observed with 
the IRS and lie in the redshift range 1.5 < z < 2.6. We 
reject the radio galaxy PKS 0529-549 (z = 2.575) as there 
is a source present in the 24 /jm MIPS map that cannot be 
distinguished from the target in the 2-d IRS data. Since the 
unknown source lies within the target extraction window, 
this makes the IRS spectrum of PKS 0529-549 unreliable. 
Our sample includes the IRS observation of 4C 23.56 pre- 
viously presented by ISevmour et alj (|2008l ). though we use 
a more rigorous spectral reduction pipeline here. We do not 
include TXJ 1908+7220 in the sample as its higher-redshift 
(2 = 3.53) results in the 9.7 /jm silicate feature and a sig- 
nificant fraction of the 7.7 /im PAH feature being shifted 
outside of the IRS spectral coverage. The properties of our 
sample are given in Table Q] and we show the 24 /im fluxes 
as a function of redshift for the whole SHzRG sample in 
Figure [T] The IRS observations are from multiple propos- 
als, with each proposal having its own particular science 
goals and we can treat the selection of our sample as be- 
ing a random selection from mid-IR 'bright' radio galaxies. 
All targets are mid-IR 'bright' in that they have 24 /im flux 
densities, S24 M m > 0.4 mjy which ensured that good qual- 
ity spectra (S/N>5 per resolution element) could be pro- 
duced in reasonable observing times. More sources from 
some of these proposals would have been observed with the 
IRS, if it was not for the depletion of the on-board cryo- 
gen. Four radio galaxies in our sample (3C 239, 4C +23.56, 
NVSS J171414+501530 and PKS 1138-26) have also been 
observed at sub-mm wavelengths with the Submillimetre 
Common-User Bolometer Array (SCUBA) on the James 
Clerk Maxwell Telescope (JCMT). Their sub-mm emission 
is discussed in 35] 



2.2 Observations 

The targets were observed with the Long-Low (LL) mod- 
ule (resolving power, R~ 60-130) of the IRS. While all tar- 
gets were observed in the first order of the LL module (LL1; 
20 /im < A < 40 /im) , only three were also observed in the sec- 
ond order (LL2; 14 /jm < A < 21 /xm). The use of the LL2 or- 
der ensured the spectra of the objects at lower redshift cov- 
ered the shorter end of our required rest-frame wavelength 
range of 6-10 /im. All but PKS 1138—26 were observed in 
staring-mode. In staring-mode, the IRS observes a target in 
two positions ('nods') in the slit of each order, resulting in 
four different 2-d spectral images if a target is observed in 
both the LL1 and LL2 orders (two images if the target is just 
observed in the LL1 order). Observations are then repeated 
in cycles for each slit and position. It was anticipated that 
the required S/N could be obtained for each spectrum by 
running observations of 20 cycles or fewer of 120 s duration 
for each target. 

PKS 1138—26 was observed in mapping-mode, where 
the target was observed at eight positions, in 1.67-arcsec 
steps, along the LL1 slit. This single astronomical observa- 
tion request (AOR) was coupled with eight other AORs to 
spectrally map the environment of PKS 1138—26 as part of 
a larger observation. The sampling of these eight AORs was 
different from the single AOR and it was intended that a 
separate, higher S/N spectrum be produc e d alo ngside the 
spectrum from the sing le AOR. lOgle et alj l|2012D describes 
these AORs in more detail and presents the higher S/N spec- 
trum of PKS 1138-26. 

For each target, the number of cycles varied depending 
on the exposure times required to reach the S/N threshold. 
The details of the observations are given in Table [2] 



2.3 Reduction pipeline 

For the seven targets observed in staring-mode, the LL mod- 
ule data obtained from the Spitzer Science Center (SSC) 
came in the form of Basic Calibrated Data (BCD) images 
produced by the S18.18 version of the IRS calibration soft- 
ware pipeline. Included in the pipeline were the low-level 
processes: ramp fitting, dark sky subtraction, linearity cor- 
rection, flat-fielding, droop correction, straylight/cross-talk 
correction, wavelength calibration and flux conversion. In 
addition to the SSC pipeline processing, bad pixel detection 
and cleaning, latent charge removal and background sub- 
traction were performed using our own reduction pipeline. 
The steps of thi s pipeline closely followed those outlined in 
iGouldind (|2010l ). The pipeline included masking damaged 
or unreliable pixels, either those previously flagged by the 
SSC calibration files or those deemed to be too variant (5 a) 
from their neighbouring pixels or their counterparts in other 
images. Masking such pixels meant they would not be con- 
sidered when combining the images from individual obser- 
vations. Images for each 'nod' in each order were combined 
using a clipped median (5er). Also performed on individ- 
ual images was the subtraction of the la tent charge that ca n 
remain on the detector between frames l|Teplitz et al.|[2007l ) . 

After the 2-d images for each 'nod' in each order were 
combined, the IRS CLEAN tool was used to replace any 
masked pixel values by interpolating the surrounding pixel 
values, effectively 'cleaning' the images. The next step was 
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Figure 1. Observed 24 /an flux density as a function of redshift for the whole SHzRG sample. The radio galaxies in our sample are 
represented by red triangles. Sources within the region enclosed by the red dashed lines correspond to SHzRG sources that are within 
our redshift range and 'bright' at 24 (im (S24uro >0.4mjy). Sources within the region enclosed by the black dashed lines correspond to 
SHzRG sources that are within our redshift range and 'faint' at 24 (im (S24(j m <0.4mjy). 



to subtract the background. For each combined image, a cor- 
responding background image was created from as many of 
the other combined images for that source as possible. An 
important factor in the number of combined images used to 
create the 'background' image was the presence of serendip- 
itous sources. A combined image would not be used for the 
background if it was deemed to have a serendipitous source 
near the target-observing position as this would result in a 
contaminated background for the target. IRSCLEAN was 
then run again, on the background-subtracted images. 

The final step in the pipeline was to perform the spec- 
tral extraction on the 2-d background-subtracted images. 
For this we used the Sp ectroscopy Modelin g Analysis Re- 
duction Tool (SMART) (|fflgdon et al.ll2004h . At the spatial 
resolution of Spttzer, our targets are point sources, so we 
used the optimal extraction option inside SMART and set 
the extraction window to be centred on the target position. 
This produced a 1-d spectrum for each nod in each order. 
Flux uncertainties were then estimated from the fluctua- 
tions along each row (i.e. wavelength) of the 2-d images. 
The extracted 1-d spectra were combined outside SMART 
and the uncertainties combined in quadrature. The spectra 
can then be compared to the 16 /jm IRS blue peak-up and 
24 ^m MIPS photometry. 

Regarding PKS 1138—26, due to the finer sampling, the 
background could not be estimated from the single AOR. 
Instead a median background i mage was created from one 
of the eight AORs obtained bv lOgle et all (|2012h . The im- 
ages for each position were combined using the standard 
Spttzer coadd pipeline and the background image was sub- 
tracted. The eight resulting images were then run through 



IRSCLEAN and the 1-d spectra were optimally extracted 
using the Spttzer IRS Custom Extraction (SPICE) software. 
The final spectrum of PKS 1138—26 is the median of the 1-d 
spectra. The individual spectra are shown in Figure 



2.4 The mean spectrum 

To obtain a measure of the average star-forming properties 
of our sample, we also create an un-weighted, arithmetic 
mean luminosity spectrum from the rest-frame spectra of the 
seven targets. So that the spectra have a constant weight- 
ing to the mean as a function of wavelength, we restrict the 
mean spectrum to only cover a wavelength range that all 
seven constituent spectra cover. This results in a coverage 
of 6.5-10 £tm. So that an accurate mean could be calculated, 
the original spectra are re-sampled. The re-sampling pro- 
cess is such that all seven spectra contribute to the mean 
in each resolution element and there is little degradation in 
the spectral resolution. 



3 MID-IR SPECTRAL DECOMPOSITION 
3.1 The model and fitting process 

To decompose the spectra into contributions from AGN and 
star formation, we perform spectral fitting on the mean 
and individual spectra in rest-frame luminosity space us- 
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ing a vari ant of the Decompllfjl IDL fitting procedure de- 
scribed in lMullanev et all (|201lf ). For the host galaxy com- 
ponent we use a li brary of temp l ates t hat comprise the 
IRS spectra of the iBrandl et al.l (|2006T l sample of local 
starburst galaxies. The se templates ha v e bee n re-reduced 
and are pre s ented in iMullanev et alj (|201ll ). From the 
IBrandl et all (|2006h sample, we exclude the galaxies that 
contain an AGN (Seyfert galaxies, low-ionisation nuclear 
emission-line region; LINERs, etc.) and accept only the pure 
starburst objects. This gives us 10 templates that differ 
in the strength of their PAH features and also their sili- 
cate absorption at 9.7 fim. The AGN component is char- 
acterised by a power-law with a spectral index, a, that 
is allowed to vary between 1.0 and 1.5. We place this re- 
striction as similar ranges for the mid-IR spectral indices of 
typ e-1 AGN have b een found (e.g. broad-line radio galax- 
ies. lOgle et alJlioOrj ; S eyfert galaxies. | Buchanan et al.l l2006 ; 
high-redshift ULIRGs,!!! onso- Herrero et aTB oOG and type- 
1 local AGN, IMullanev et alj|201lh and the mid-IR slopes 
(5 /im < A < 38 fim) of ty pe-1 and type-2 objec ts have been 
shown to be comparable (Ne nkova et al.ll2008al ). To account 
for the higher inclinations of our type-2 sources, we make 
an assumption that the dusty environments of high-redshift 
radio-l oud AGN can be modelled by a Galactic extinction 
curve |Saiina et alj|2007al : lLacv et all |200?1 ; ISevmour et all 
20081) which w e apply to the power-law slope. We choose the 
Draind <|2003h extinction curve that uses Rv = 3.1, the typi- 



cal value for the Milky Way Galaxy, where Rv is the ratio of 
total extinction, A(V), to selective extinction, A(B)-A(V), 
in the V band (i.e. R v = A(V)/(A(B)-A(V))). Although this 
is a simplistic approach as we assume that the obscuring 
medium acts as a simple foreground screen, it will suffice for 
fitting AGN continua and silicate absorption features in our 
spectra. We d o not apply extincti on to the host galaxy com- 
ponent as the lBrandl et al.l (|2006T ) spectra are empirical and 
have a certain level of obscuration already, evident through 
their silicate absorption features. 

The five free parameters for our model are: (i) the host 
galaxy template, (ii) the host galaxy normalisation, (iii) the 
power-law normalisation, (iv) power-law spectral index, a 
and (v) the hydrogen column density, JVh- The expression 
for the modelled mid-IR spectrum, L„, in WHz -1 is thus: 



L v (\) = ki X a e 



-N H a(A) 



k 2 fu(X) 



(1) 



Where f v (A) is the flux density, in Jy, of the host galaxy tem- 
plate, ki and are the normalisation factors, and cr(A) is 
the extinction cross-section that, with the hydrogen column 
density, gives the optical depth, t(A); i.e. t(A) = A^h o"(A). 

For each host galaxy template, x 2 minimisation is per- 
formed on the spectra, simultaneously solving for fci, k-z, a 
and Nh- In Figures [2] and [3] we show the resulting best- 
fitting components for the individual spectra and mean 
spectrum, respectively. To determine the confidence in the 
X 2 minimisation, the free parameters are allowed to devi- 
ate from their best-fitting values until the condition x 2 ~ 
Xmin — A-X 2 = a 2 is met; letting a = 1 or 3 to find 1 a and 3 a 
uncertainties, respectively. The extent of the free parameter 



1 Our variant is that the procedure only involved IRS spec- 
troscopy; i.e., no broad-band IR photometry was fitted. 

2 https:/ /sites. google.com/sitc/dccompir 



deviations from the best-fitting values represents their un- 
certainty limits. If Ax 2 < 9 is returned when the host galaxy 
component is not included, then at the 3 a detection level, 
the fit is consistent with a scenario in which there is no star 
formation. 



3.2 Deriving host galaxy and AGN properties 

For the templates that fit best, we use the host galaxy nor- 
malisation fa c tors t o convert the 7.7 fim PAH fluxes given in 
IBrandl et al.l (|2006T ) into 7.7 fim PAH luminosities, L7.7PAH, 
for our sample. The range of acceptable normalisation fac- 
tors gives the range in PAH luminosity and thus, the un- 
certainty on the best-fit value. The PAH luminosities are 
then converted to total IR (8-1000 fim) luminosities, im,, (if 
assumed to be entirely due to star formation) through the 
conversion in P08 which is given as: 

logL 7 .7PAH = (-1.5 ± 0.1) + (0.95 ± 0.01) log Lir (2) 

This conversion was determined empirically by P08, through 
the PAH and sub-mm luminosities of a combined sample of 
local starburst galaxies and high-redshift, star formation- 
dominated sub-mm galaxies (SMGs). These SMGs lie at 
similar reds hifts to our sou ces. Lir, is then converted to 
SFR via the lKennicuttl (|l998T ) relatio n below, which a ssumes 
a starburst period < 100 Myr and a ISalpeterl (| 19551 ) initial 
mass function (IMF): 

l > Lir 



1 M yr- 



5.8 x 10 9 L e 



Furthermore, the sample of P08 exhibit a 1 a scatter of a fac- 
tor of 2 about the derived L7.7PAH — Lir conversion which 
translates to a range of SFR estimates from each derived 
PAH luminosity. The effect of this choice of conversion fac- 
tor is investigated in £14.11 where we consider the conversion 
derived by MD09, in which a different IRS sample of SMGs 
was used. 

From the AGN component, we determine monochro- 
matic luminosities at rest-frame 6 fim, vL v (Qfim), where 
the continuum is less contaminate d from spectral fe atures 
than it is at longer wavelengths (|Lutz et al. 2004). The 
spectrum for PKS 1138—26 and also the mean spectrum 
(due to reduced coverage) do not cover rest-frame 6 fim and 
so when estimating their 6 fim luminosities we assume the 
AGN emission at this wavelength is still characterised by 
the measured power-law and this component is simply ex- 
trapolated. In both cases, the extrapolation is < 0.5 fim. We 
also quantify the strength of silicate absorption at 9.7 fim 
i.e. the 9.7 fim apparent optical depth, by taking the natu- 
ral log of the ratio of continuum emission, / CO nt , to observed 
emission, / G bs, i.e. T9.7 M m = ln(f cont /fob B ) = AHO-g.7 Mm ; 
09.7/jm ~3.9 x 10 _23 cm 2 . Note that we assume a clumpy 
rather than homogeneous distribution for the obscuring 
medium (see ^SJ so our definition of optical depth is no t 
equivalent to the traditional one (e.g. iNenkova et al]|2008al ) . 
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Figure 2. The mid-IR spectra of our sample along with their decomposition into AGN and host galaxy components. For all spectra, the 
data and the 1 a uncertainties are represented by a histogram and shaded region, respectively, while the overall best-fitting is denoted by a 
solid red line, the host galaxy component is denoted by a dotted green line, while the power-law and extinguished power-law components 
are denoted by dashed blue and dot-dashed orange lines, respectively. Broadband MIPS/IRS photometry, though not included in the 
fitting, is also plotted (filled purple circles). The location of potential PAH features are indicated by arrows. The [SIV] emission line 
at 10.5 /Ltm is visible in the spectrum of 3C 239 and PKS 1138-26 while the [Arlll] line at 9yum is visible for PMN J0408-2418. The 
location of the [S IV] and [Ar III] lines in the three spectra are indicated by vertical dashed and dot-dashed lines, respectively. 
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Figure 3. The mean spectrum of the seven radio galaxies in the sample, with the same fitted components as those in Figure [2] The 
location of potential PAH features are indicated by arrows. [Arlll] is detected at 9 fim and its location is indicated by a dot-dashed line. 
The location of [S IV] lies outside the spectral coverage. 



4 RESULTS 

4.1 Spectral and fitting results 

The individual spectra are presented in Figure [2] Note that 
the continua are consistent with the 16 /jm IRS blue peak-up 
and 24 fim MIPS photometry. The mean S/N per resolution 
element for the whole sample ranges from ~ 3-23. The mean 
spectrum has a mean S/N per resolution element of 37 and 
is shown in Figure [3] 

We find that all sources are AGN-dominated in the mid- 
IR, which is in agreement with the pr evious IR SEP mod- 
elling of the whole SHzRG sample (e.g. lSevmour et alj|2007t 
DB10). However, 6/7 also show a significant host galaxy con- 
tribution. For the other source (NVSS J180556+633313), a 
fit that includes no host galaxy component differs from the 
best-fitting model by A\ 2 < 1, so we do not significantly 
detect star formation in this source. For the radio galaxies 
whose spectra have a mean S/N per resolution element < 10 
(NVSS J180556+633313, NVSS J171414+501530 and 4C - 
00.54), more than one host galaxy template can be fitted 
with A\ 2 < 1. In this case, similar 7.7 /im PAH luminosities 
are derived from the acceptable templates and so although 
the PAH luminosities have been derived from a model-fitting 
process, they are at least insensitive to the choice of host 
galaxy template. 

At the extreme, we calculate very high SFRs of 
~lOOOM yr _1 (for 4C 23.56, PMN J0408-2418 and PKS 



1138-26). We find an SFR of 620 t ™ M yr* 1 for the mean 
spectrum and a mean SFR of 610 ± 180 Mq yr _1 for the in- 
dividual sources. The error of the mean SFR is derived from 
the scatter of the individual SFRs and does not take into 
account the PAH to SFR conversion uncertainty. We exam- 
ine how the radio galaxy with the highest equivalent width 
7.7 /xm PAH feature in its spectrum, PMN J0408-2418, con- 
tributes to the mean SFR by creating a mean spectrum 
without including the source. We derive a mean SFR of 
510 ^260 Mq yr _1 suggesting that star formation is abun- 
dant in our sample. Except for NVSS J180556+633313, the 
total IR luminosities (inferred from 7.7 /im PAH luminosi- 
ties) alone qualify these objects as ultra luminous infra-red 
galaxies (ULIRGs; 12 ^ log(LiR,/LQ) < 13). Assuming the 
IR-radio correlation for star-forming galaxies hol ds at high- 
redshift l|Sevmour et alj|2009l ; Irluvnh et alj|20ich , we calcu- 
late star- form i ng ra dio luminosities from the derived SFRs 
using the Bell (2003) conversion factor. The star-forming ra- 
dio luminosities are ~ 5 orders of magnitude lower than the 
measured radio luminosities and so the contribution of star 
formation to the radio emission in these HzRGs is negligible 
as expected. We next investigate the effect our choice of the 
P08 -L7.7PAH — Lm conversion has on our results. Applying 
the MD09 conversion (see £|3.2[1 to the mean spectrum we 
find a very similar SFR of 710 ± 40 M yr -1 , well within the 
1 a uncertainties using the P08 conversion. No 1 cr scatter in 
the relation is given by MD09 and so the uncertainty given 



© 2012 RAS, MNRAS 000, [THM] 



8 J. I. Rawlings et al. 



here reflects only the uncertainty in the PAH luminosity. 
The similarity between the SFRs from the two conversion 
factors is an indication that our SFRs are not sensitive to the 
choice of conversion factor from PAH luminosity. To derive 
their L7.7PAH —Lm relations, both P08 and MD09 used sim- 
ilar methods and objects to us, in that their mid-IR spectra 
were decomposed into AGN and host galaxy components. 
Therefore, it is appropriate to apply these conversions for 
our sources. 

Using the stellar masses from Table [T] and the SFRs 
from Table [3] we now calculate the specific star formation 
rate (sSFR) of our sample. The sSFR is a measure of a 
galaxy's rate of growth with respect to its stellar mass. We 
do this for the two sources that have stellar mass estimates 
and we find lower limits for the other sources that only have 
upper limits for their stellar mass (due to unknown con- 
tributions from the AGN to the //-band luminosities, from 
which the masses are derived) . The exception to this is NVSS 
J180556+633313, for which we only have upper limits for the 
SFR and stellar mass and therefore are unable to calculate 
its sSFR. It is noteworthy that different IMFs have been 
used to c alcuat e the stellar masses and SFRs (see H3.2|) . The 
iKroupal l|200lf ) IMF underestimates the total stellar mass 
compared to the lSalpete r (1955) IMF, for a galaxy of a given 
luminosity, by a factor of ~ 1.5 (Karczewski et al. 2012; in 
prep.). We account for this factor when calculating the sS- 
FRs for our sample. We derive sSFRs of 0.31 ± ° l 2 7 Gyr -1 
and 3.12 tr.sg Gyr" 1 for 3C 239 and PMN J0408-2418, re- 
spectively, while three of the other four galaxies have sSFRs 
> lGyr" 1 . 

From the AGN component, we find very high mid-IR lu- 
minosities, with vL v (6 (im) « 10 46 ergs -1 for the mean spec- 
trum and vL v (6/im) between 10 45 ' 2 and 10 46 ' 5 ergs -1 for 
the individual objects. We also find that similar AGN lumi- 
nosities were derived when each template was fitted, show- 
ing that the AGN component is not sensitive to the choice of 
template. As for the obscuration of the AGN component, we 
find low to moderate 9.7 /im optical depths (T9.7 Mm < 0.8) 
for the sample. When investigating the degeneracy of a and 
extinction, we find little covariance of the 1 a confidence 
regions in the parameter space and so we observe little de- 
generacy. 

Although wc find many of these radio galaxies have sig- 
nificant star formation, one important point to note is that 
the sample could be biased towards objects with significant 
star formation because, at z ~ 2, the 7.7 /im PAH feature 
shifts into the 24 /im MIPS band. While these radio galax- 
ies were not 24 /im selected, they are all bright at 24 /im 
(S24 M m >0.4mJy) as this was a requirement for IRS obser- 
vations. This means a substantial fraction of the 24 /im emis- 
sion could be from star formation. In Table O we consider 
the contribution of the host galaxy component to the total 
flux in the 24 /im MIPS passband for our sample. For approx- 
imately half the sample, these contributions are significant 
(> 20 %) but only one source would potentially drop out of 
the sample if the star- forming contributions were subtracted 
from the 24 /im flux densities. Therefore, these sources are 
mainly selected on their AGN power but the star-forming 
bias may still exist. 
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Figure 4. 7.7/tm PAH luminosity as a function of mid-IR AGN 
luminosity for the sample. The mean spectrum is represented by 
a star. 

4.2 Average star formation rate of the HzRG 
population 

As mentioned previously (in j)2.1|) . the original SHzRG sam- 
ple is representative of the HzRG population. In addition, 
amongst the sources that are within our redshift range and 
bright at 24 /im (S24 Mm > 0.4 mjy) i.e. the sources within the 
region enclosed by the red dashed lines in Figure [TJ we can 
treat those selected for IRS observations as a random selec- 
tion. As a result, even though we have estimated the SFR of 
only seven HzRGs, we can estimate the total mean SFR for 
the overall HzRG population in our redshift range (1.5 < 
z < 2.6). This is possible with a few assumptions. We can 
assume the mean SFR (( SFR > = 610 ± 180 M yr -1 ) of our 
sample is equal to the mean SFR of the sources within the 
region enclosed by the red dashed lines in Figure [T] We can 
also make two assumptions on the SFR of the sources in our 
redshift that are 'faint' at 24 /im (S24,jm < 0.4 mjy) i.e. the 
sources within the region enclosed by th e black dashed lines 
in Figure [1] As in lSevmour et all (|201ll ). firstly, to obtain a 
lower limit to the total mean SFR, we assume that the mid- 
IR faint sources have SFRs of zero. The total mean SFR 
then just scales with the mean SFR of the mid-IR-bright 
sources by the fraction of mid-IR bright sources (i.e. to- 
tal mean SFR = 24/30 x ( SFR ) = 490 M yr" 1 ). Secondly, 
to obtain an upper limit, we assume that the mid-IR faint 
sources have SFRs equal to the mean SFR of the mid-IR 
bright sources. With this approach, we constrain the total 
mean SFR of the HzRG population in our redshift range to 
be 490-610 Mq yr -1 . Taking into account the standard error 
of the mean of ± 30% gives the range 340-790 M Q yr -1 . 

4.3 Relating AGN to host galaxy parameters 

We now investigate the relationship between AGN and host 
galaxies, conducting Spearman's rank statistics to identify 
possible correlations for the physical parameters of our sam- 
ple. More specifically, we determine the Spearman's rank co- 
efficient, p, and the significance of the deviation of p from 
zero, a. Values very close to unity for p and a would indicate 
a significant correlation exists. 

In Figure [4] we plot 7.7 /im PAH luminosity against 
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Figu re 5. SFR as a function of 500 MHz luminosity for our sources (filled black circle s), the high -redshift radio-loud AGN from 
ISaiina et all d2007bl , blue asterisks) and TXJ 1908+7220 (2 = 3.53, orange triangle, [Seymour et al.ll2008h. Also plotted is t he mean SFR 
for the 30 SHzRGs objects in our redshift range (grey region) and the p arameter space cover ed by the lSevmour et al.l ll201ll ) high-redshift 
(1.2 < z < 3) sample of radio-loud AGN (green hashed region). For the lSaiina et al.1 j2007bl ) sources, 500 MHz luminosities are calculated 
from 1.4 GHz lumino sities using their radio spectral indicies. For TXJ 1908+7220, we use the 500 MHz luminosity from DB10. For the 
ISevmour et al.l i201ll ) source s, 500 MHz luminosities are calculated from 1.4 GHz luminosities by assuming a radio spectral index of -0.75. 
For the TSa!inaet**aT[ l|2007bl ) sources and TXJ 1908+7220, SFRs are calculated from 7.7 /im PAH luminosities using the same method 
as our sample. The mean spectrum is represented by a star. 



rest-frame 6 /im AGN luminosity for all our sources and 
the mean spectrum. The scales for both parameters vary 
by more than an order of magnitude, showing a large di- 
versity in the limited sample and we find only a marginal 
correlation (p = 0.67, a — 0.90). Next we plot SFR versus 
500 MHz luminosity for our sample and the mean spec- 
trum in Figure [5] along with results from the literature. 
We also show the mean SFR for the 30 SHzRGs objects in 
our redshi ft range (grey r egion), the high-redshift radio- loud 
AGN from ISaiina et al.l (|2007bl . blue asterisks), the higher- 
redshi ft radio galaxy TXJ 1 908+7220 {z = 3.53, orange tri- 
angle, ISevmour et al.l l200ct ) a nd the range in mean SFR 
and 500 MHz luminosity for the lSevmour et al.1 (|201ll ) high- 
redshift (1.2 < z < 3) sample of radio-loud AGN (green 
hashed region). It is visible that mid-IR bright radio- loud 
AGN over a few orders of magnitude in radio luminosity 
can have significant star formation. 

Radio core dominance, R, defined as the ratio of core 
and extended radio flux at rest-frame 20 GHz, is a measure 
of the extent to which the central component of the radio 
source contributes to the total radio emission. It can there- 
fore be used as a proxy for the orientation of the radio jets 
|Kapahi fc Saikialll982l . Drouart et al. 2012; in prep.). Using 
the radio core dominances from Table [TJ we next consider 



if jet orientation shows any relationship with the mid-IR 
obscuration of these AGN. Due to anisotropic, relativistic 
beaming effects, the cores of the radio sources will appear 
brighter, relative to their jets, if the jets are orientated close 
to the line of sight. In Figure[S]we plot 9.7 fim optical depth 
versus the radio core dominance for our sample and observe 
no correlation. A correlation may have been anticipated as 
AGN mid-IR obscuration and jet orientation are expected 
to correlate if their axes are perpendicularly alligned. The 
likely reason for observing no correlation is discussed in Sj5] 

4.4 Other spectral features 

The rest-frame 10.5 /im [S IV] atomic emission line is de- 
tected in the spectra of 3C 239 and PKS 1138-26, 
with equivalent widths (EWs) of 0.07 + 0.01 /im and 
0.04 ± 0.01 /im, respectively. For PMN J0408-2418, the rest- 
frame 9 /im [At III] emission line is detected with an EW of 
0.06 ± 0.01 fim. We also detect [Arlll] in the mean spec- 
trum (EW = 0.015 + 0.003 ^m) but the location of [SIV] 
lies outside the spectral coverage. These EWs are larger 
than those found for a sample of loc al Seyfert galaxies 
(EW [S iv] <0.03 /im, EW[a,ttt] < 0.01 um. lffimig et al.ll2008h 
and are comparable to those for other powerful type-2 AGN 
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Figure 6. 9.7 (im optical depth as a function of radio core dom- 
inance. The mean spectrum is represented by a star. 



i|Shan fc CherJl2012h . No H2 emission is detected in any 
of our spectra. However, for PKS 1138—26, H2 is detected 
in the higher S/N IRS spectrum presented in lOgle et all 
i|2012h . Therefore, the non-detections in our spectra may be 
partly due to low S/N. 



5 DISCUSSION 

Through the detected PAH emission in their mid-IR spec- 
tra, we observe that some HzRGs are forming stars at a very 
high rate. We find a high mean SFR of 340-790 M Q yr" 1 for 
the overall HzRG population in the redshift range 1.5 < 
z < 2.6. However, not all HzRGs at this epoch have 

the case for NVSS 
High SFRs have 



such significant star formation, as is 



J180556+633313 (SFR = 20 + f M yr' 
also been detected in other high redshift radi o-loud AGN, 
for exa mple, the AGN-dominated ULIRGs of ISaiina et al.l 
(|2007bl . SFR«lOOOM yr -1 ). However, these sources were 
selected to be bright in the mid-IR (S24 M m > 0.9 mjy) before 
they were discovered to be radio-loud, implying that the 
sample might have been biased towards objects with signif- 
icant star formation. As show n in Figures [4] (right panel) 
and \5\ the lSaiina et~aH |2007bl ) sample have lower radio lu- 
minosities than our sources. Therefore, PAH emission exists 
in mid-IR bright radio-loud AGN over ~ 4 orders of magni- 
tude in radio luminosity. In Figure [5] we also plot the less - 
luminous high-redshift AGN from ISevmour et al.l (|201l[ ). 
which have a mean SFR of 80-580 M yr" 1 , comparable to 
or less than the mean SFR for the HzRG population in our 
redshift range. For our sample, we have not found a signif- 
icant correlation between star formation and AGN power, 
for both mid-IR and radio luminosities, of po werful radio- 
loud AG N at high redshift. When i ncluding the lSaiina et al.l 
(2007b]) and ISevmour et alj l l201ll) samples, the scatt er in- 
creases. This is consistent with iMullanev et al.l (|201 If ) who 
also examined SFR as a function of AGN power. However, 
other studies have been able to show a connecti on between 
star f ormation and centra l black hole accretion l)Shao et al.l 
l2010l : lBonfield et al.l[201lh . 

To further investigate whether the use of the 
£<7.7pah — Lm conversion from P08 was appropriate 



Figure 7. Ssso fim/i7.7PAH verses redshift for the sources in our 
sample with sub-mm observations (black circle/upper limits) and 
the sub-mm galaxies from P08 (red triangles) and MD09 (blue 
squares/lower limits). The lower limits for some of the MD09 
sources are due to greater contributions from noise in the 7.7 /an 
PAH region. 



for our sources, we can compare the PAH and sub-mm 
emission from our sources to the emission from the IRS- 
observed SMGs that were used to (partly) derive the 
conversion. As mentioned previously, four radio galaxies 
in our sample have sub-mm observations; three have non- 
detections (3C 239, S 8 5o M m= 0.00 ± 1.00 mjy; 4C +23.56, 
S850 M m = 1-68 ± 0.98 mjy and NVSS J171414+501530, 
Sgso^m = 0.99 ± 1.10 mjy; corrected for radio sy nchrotron 
emission contamination; lArchibald et al.l [2001) and the 
other is detected (PKS 1138-26, Sgso^m = 6.7±2.4mJy; 
IStevens et all [20031. For these sources and the SMGs from 
P08 and MD09 that had full spectral coverage of the 7.7 /im 
PAH feature (12/13 and 21/24 objects, respectively), 
we plot the ratio of 850 /im flux density to 7.7 /im PAH 
luminosity as a function of redshift in Figure Due to 
the strong negative ^-correction at 850 fim for thermal 
dust emission, the 850 fim flux density can be used as a 
crude proxy for far-IR luminosity at z > 1. We retrieve d 
850 /j.m flux den si ties fo r the SMGs fromlPope et all |2006h . 
IChapman et all (|200ll 1 and IChapman et alj l|2005h . while 
redshifts were either derived from the IRS data by P08 
and M D09 or optical s pectro scopic redshifts were retrieved 
from IChapman etafl (|2005h . We find our radio galaxies 
have S850Mm — i7.7PAH ratios (upper limits for the sub-mm 
non-detections) within the ranges of both samples of SMGs 
and so this indicates that it was reasonable to apply the 
i7.7PAH — Lm P08 conversion to our sources. 

We turn now to the 9.7 silicate features in the IRS 
spectra. We begin by examining whether the host galaxy 
could contribute significantly to the 9.7 ^m silicate feature 
by considering the emission at shorter wavelengths where the 
effects of extinction are greater than in the mid-IR. Near-IR 
(rest-frame optical) spectra exist for five sources in our sam - 
ple; namely 4C 23. 56 and 4C-00.5 4 dHumphrev et alj2008h. 
PKS 1138-26 l|Nesvadba et~all 120061; iHumphrev et all 
l200ct l. NVSS J171 414+501530 jMotohara et al.ll200lT ) and 
PMN J0408-2418 (|McCarthv et al.lll992T ). Only one source 
(PKS 1138—26) has a significant broad component to Ha, 
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from which iNesvadba et al.l (|2006l ) estimated the total ex- 
tinction in the V-band, Ay. They find Av = 8.7, and this 
source has a very weak silicate feature {Tg.r^m < 0.07 at 
3<r). Whether the host galaxy contributes to the feature 
or not, any silicate absorption from the torus must be weak 
(i"9.7/jm < 0.8) because absorption from the host galaxy can 
only increase the dep ths of the feature s . From the narrow- 
line Ha / H/3 ratios, iHumphrev et all (120081 ) placed con- 
straints on the Av of the extended emission line region 
(EELR) for 4C 23.56, PKS 1138 -26 and 4C-00.54 These 
regions are external to the torus llFosburv et al.ll 19821 . [l98i 
iBaum et al.l Il988l ; iTadhunter et al.l Il989h so any obscura- 
tion is likely to be due to the host galaxy. Both 4C 23.56 
and PKS 1138—26 have narrow-line Ha /H/3 ratios con- 
sistent with case B photoionisation with no extinction. A 
greater Ha / H/3 ratio was derived for 4C— 00.54, suggesting 
non-neglig ible extinct i on fro m the host galaxy. From the ra- 
tios, |Hu^hre^elaL| (120081 ) derived Ay < 1.2 and < 0.7 for 
4C 23.56 and PKS 1138-26 respectively, and Ay > 1.2 for 
4C— 00.54. We translate Ay to silicate absorption strength of 
the EELR, T9.7 Mm (EELR), using the same extinction curve 
used in the fitting and find r 9 . 7 M m (EELR) < 0.09, < 0.05 
and > 0.07 for 4C 23.56, PKS 1138-26 and 4C-00.54 re- 
spectively. These values are consistent with the absorption 
strengths from the fitting of to,.t ^ =0.12 ±0.06, < 0.07 and 
< 0.51 (3<r upper limits). Regarding the near-IR spectra for 
NVSS J171414+501530 and PMN J0408-2418, narrow Ha 
emission is detected but there is no coverage of the H/3 line 
and so we cannot constrain the host galaxy extinction for 
these sources. Considering the sample overall, we cannot rule 
out a contribution from the host galaxy extinction to the 
silicate absorption features, but neither can we identify any 
silicate absorption features for which we are confident that 
the host galaxy contributes. 

To characterise the tori in our sources, we next con- 
sider the attempts that have been made to model the 
torus geometry using the 9.7 /im sili cate feature. It ha s 
been shown t hat both homogeneous JPier fc Krolik 1992) 
and clumpy (|Levenson et al.l 20071 : iNenkova et al. 2008b) 
torus models can synthesise the silicate feature observed 
for a number of type-2 A GN |Hao et al.ll2007l ; ISpoon et all 
120071 : iNikutta et atll2009l ). However, a homogeneous geom- 
etry is unphysical for the torus due to thermal spallation 
jKrolik fc Begelmarj|l988l ) and although such a geometry is 
a reasonable approximation as long as the distance between 
the clouds is small in compari son to the scale of the whole 
structure (|Pier &: Krolikl Il992|). we consider only a clumpy 
distribution from here on in. iNenkova et~aH (2008a) showed 
that for a clumpy torus, a silicate feature in emission could 
be produced, particularly if there are very few clouds along 
any radial equatorial ray. Such emission from the torus could 
potentially be masked by stronger silicate in absorption from 
the host galaxy. We cannot rule out such a situation from the 
data but as we have no evidence to support this, we make 
the caveat that in the following discussion we attribute the 
silicate features (or lack of) entirely to the torus component , 
in line with other similar studies (e.g. Levenson et al.ll2007i : 
ISpoon et all 120071 ; ISaiina et all l2007al: IClearv et all l2007l) . 
Considering the model from Nenkova et al.l ( 2008a ) , we con- 
strain some of the p hysical parameters of t he torus for our 
AGN. For instance, INenkova et alj (|2008al ) showed that a 
soft-edged distribution better matches observations of the 



nuclear emission of a sample of Seyfert galaxies than a sharp- 
edged geometry does (see their Figure 1 for a schematic di- 
agram of the two toroidal distributions). Furthermore, their 
models were able to reproduce the 9.7 /im silicate feature 
properties of Seyfert galaxies when the average number of 
clouds along any radial equatorial ray is ~5-15, the angu- 
lar width of the torus is ~ 15-45 and the actual optical 
depth of the individual clouds is ~ 30-100. Considering PMN 
J0408— 2418; the one HzRG for which we calculate a moder- 
ate 9.7 /im optical depth (rg.7 ^ = 0.70 1 o'r?), f° r the torus 
model to reproduce this absorption strength, an average of 
~ 10-15 clouds along any radial equatorial ray are required, 
with each cloud having an optical depth of ~ 10-80. An angu- 
lar width of ~ 30-45 is also needed. This particular model 
has a radial distribution profile of r -2 and an outer-to-inner 
radius ratio of 30. A model with an angular width as low 
as 15 ° can reproduce the absorption strength but for this 
a cloud optical depth of ~ 20-50 is required. For the other 
HzRGs that have weaker 9.7 /jm optical depths, the silicate 
absorption strengths can be reproduced by a model that has 
an angular width of ~ 15-45 ° and also a radial distribution 
profile of r~ 2 and an outer-to-inner radius ratio of 30. How- 
ever, we cannot differentiate a model that has a lower num- 
ber of clouds, ~ 5, and a moderate cloud optical depth, ~ 40, 
with a model that has a higher number of clouds, ~ 10-15, 
and either a much lower or much higher cloud optical depth, 
~ 10 or ~ 100. This is because both the number of clouds 
and the cloud optical depth are the major driving forces 
behind the strength of silicate absorption (see Figure 6 in 
INenkova et al.ll2008al ). implying a higher level of degener- 
acy in these fitting parameters when modelling the torus of 
AGN with weaker optical depths. 

To investigate the connection between AGN mid-IR ob- 
sc uration and o ri entat ion, we compare our results to those 
of IClearv et Zl (|2007f ). For an orientation-unbiased sam- 
ple of lower-redshift (z < 1) radio galaxies, they find an 
anti-correlation between 9.7 /im optical depth and radio core 
dominance which implies there is a coupling of the axes of 
the radio jets and the mid-IR obscuration. Although we ob- 
serve no correlation between these parameters in Figure [6l 
for the same range in core dominance, our sources cover a 
similar range in optical depth. The reason we do not ob- 
serve a correlation is likely to be because the smaller sam- 
ple of our study means we probe a smaller region in pa- 
rameter space which makes it more difficult to observe the 
correlation. Comparing to other studies, our 9.7 /tm opti- 
cal depths are s i milar to those for radio-quiet type-2 AGN 
(|Hao et al.ll2007l;ISpoon et al.ll2007h and other low er-redshift 
radio galaxies (|Haas et al.ll2005l ; IOgle et al.ll2006l ). This sug- 
gests that the tori in these AGN have similar properties to 
the tori in our sample. However, on average, higher levels 
of obscuration (79.7 Mm > 1) have been found fo r other sarn- 
ies, su ch as the radio- loud AGN/ULIRGs of ISaiina et al.l 
2007b|) and local Compton -thick (/V H > 1.5 x 10 24 c m~ 2 ) 
AGN ((Goulding et al.|[2012l ). For lSaiina et"aH |2007bl ), this 
is likely due to their sources being selected to have large mid- 
IR-optical colours. The authors from both studies attribute 
their larger optical depths to host-galaxy contribution to 
the obscuration of the AGN by way of acting as a cold dust 



We next examine these radio galaxies in the con- 
text of galaxy formation and evolution. The sSFR of 
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0.31 1 p'^ 2 Gyr 1 we calculate for 3C 239 is co nsistent with 
local star-forming galaxies fsSFRss 0.25 Gyr -1 . lElbaz et al.l 
2011), however, most of the other sources have much higher 
sSFRs and so they are building their mass at a faster 
rate. Assuming a constant SFR, most of the radio galax- 
ies have taken < 1 Gyr to form their stars and are likely 
to have already formed the bulk of their mass. This is be- 
cause galaxies with masses of ~ 10 12 Mq reside in the tail 
of the stellar mass func t ion u p to at least z = 4, as shown 
in lPerez-Gonzalez et al.l (|2008l ). If the star formation in our 
sources continues at the rate we observe, it would take them 
~4Gyrs to reach the ste llar mass of the gian t elliptical 
galaxy M87 (~ 10 12 ' 5 M , IWu fc Tremaindl2006l ). If we as- 
sume the mass of M87 as an approximate limit to a galaxy's 
stellar mass, it is unlikely that our radio galaxies will main- 
tain their current SFRs to the present day. A possible sce- 
nario for quenching t he star formation is t he evolutionary 
sequence proposed by I Sanders et all (| 19881 ). This sequence 
is as follows: (i) through a major-merger event, there is a 
phase of extreme star formation, (ii) once a supply of gas 
has reached the central region, an obscured AGN arises, (iii) 
as the SFR declines, the AGN starts to dominate, it sheds its 
obscuring dust and a quasi-stellar object (QSO) is formed, 
(iv) once the QSO runs out of fuel, this culmin ates into a qui- 
escenc e elliptical galaxy. The first stage of the lSanders et al.l 
(1988) sequence is consistent with the proposition that pow- 
erfu l radio-loud AGN ar e a consequence of gas -rich merg- 
ers (jHeckman et al.lll986l ; lHardcastle et al.ll2007f l. After the 
QSO phase, it is now thought that the AGN inhibits further 
star formation through 'radio-mode' feedback, in which the 
relativistic jets inject enough energy into the intergalactic 
medium to suppress the cooling of gas onto t he AGN host 
galaxy and the eventual production of stars (I Bower et al.l 
120061 ; ICroton et alJlioO^ ). The inclusion of SMBH activity 
in galaxy formation models therefore helps to explain why 
most local massive elliptical galaxies predo minantly con- 
tain o ld and red stars. In the context of the lSanders et all 
( 1988) evolutionary sequence, the properties we measure for 
our sample of radio galaxies (high AGN luminosities and 
SFRs) could be the result of a merging event, however we 
no w observe them at a p oint early on in the third stage of 
the I Sanders et al.l fl988) evolutionary sequence, where the 
obscured AGN dominates the system but star formation 
is still significant. The radio galaxy 3C 239 (and perhaps 
NVSS J180556+633313), with a relatively high stellar mass 
and low SFR, appears to be further along this evolutionary 
track than the rest of the sample, where a QSO is closer to 
forming. Furthermore, the short timesc ales su ggested f or the 
radio -loud phase of AGN (< 10 7 years. iBlundell k, Rawlings] 
GUI) in combination with the high stellar masses of our 
HzRGs, implies that they must have been forming their stars 
long before the onset of the radio emission. This suggests 
that the powerful radio-loud phase occurs near the end of the 
formation of the host gala xy, perhaps du r ing th e third stage 
of evolution proposed by ISanders et all (| 19881 ). but before 
the point at which host galaxy activity ceases completely, 
allowing the two phases to be observed simultaneously. 



6 CONCLUSION 

We have separated the mid-IR spectra of seven power- 
ful high-redshift radio galaxies (1.5 < z < 2.6) into 
AGN and host galaxy components, observing PAH emission 
and silicate absorption for several sources. These HzRGs 
are powerful systems with high mid-IR AGN luminosities 
(yh v (6 /im) « 10 46 ergs _1 for the mean spectrum). 

We find a wide range of SFRs, ranging from 
~ 1000 M© yr _1 to no significantly detected star formation 
and we constrain the mean SFR for the overall HzRG 
population in our redshift range (1.5 < z < 2.6) to be 340- 
79OM0yr~ 1 . Thus, at the epoch we probe, some HzRGs 
are forming stars at a high rate. Comparing with the results 
from other studies, we observe that radio-loud AGN over a 
few orders of magnitude in radio luminosity have comparable 
SFRs. Overall we find no evidence for a significant correla- 
tion between radio emission and star formation in powerful 
radio- loud AGN at high redshift. 

Attributing the silicate absorption feature entirely to 
the AGN torus, the relatively small absorption features 
(T9.7 Mm < 0.8) we derive for these objects is an indication 
that their tori have a clumpy distribution and contain rela- 
tively small number of clouds, ~ 5-15, through the equatorial 
plane. These optical depths agree with those for some radio- 
quiet type-2 AGN and lower-redshift radio galaxies, imply- 
ing the tori have similar properties. However, the greater 
silicate absorption features observed for radio-loud AGN- 
dominated ULIRGs and local Compton-thick AGN suggests 
the host galaxy can act as a cold dust screen for some 
sources. 

The sSFRs we calculate range from ~0.3 to 3 Gyr -1 , 
although lower limits of > 1 Gyr -1 are calculated for most 
of the sources in the sample. Assuming a constant SFR, this 
suggests that most of the sources have taken < 1 Gyr to form 
their stars. They are also likely to have already formed the 
bulk of their stellar mass and so we expect their SFRs to 
decline in the future. 

To consider how these radio galaxies formed and 
evolved, they likely represent a stage where the obscured 
AGN has developed enough to dominate the system but the 
declining star formation is still significant. 

Since the timescale of the radio emission is short and 
these radio galaxies are all massive, they must have started 
to form their stars a long time before the radio-loud phase 
commenced. This is an indication that the powerful radio- 
loud phase occurs near the end of the host galaxy activity 
but during a stage where significant star formation is still 
taking place. 
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Table 1. The properties of the seven HzRGs in the sample. Columns are as follows: (1) HzRG name; (2) name as given by the SHzRG project; (3) and (4) J2000 co-ordinates; (5) 
spectroscopic redshift; (6) rest-frame 500 MHz luminosity, with assumed uncertainty of 20% ; (7), (8) and (9) observed 16 fim, 24 /im and 850 (im flux densities, respectively; (10) stellar 
mass; (11) angular size of radio source, with assumed uncertainty of 10%; (12) radio core dominance, R, defined as the ratio of core and extended radio flux at rest-frame 20 GHz, with 
assumed uncertainty of 50%. 
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[WHz- 1 ] 




[mjy] 


[mjy] 


[mjy] 
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(11) 
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1 AC O O 1 C\ 
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< 11.59 


53.0 


0.171 


NVSS J 180556+633313 
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NVSS J171414+501530 
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+50:15:30 
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PMN J0408-2418 
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-24:18:16 
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2.427 


29.03 





64 + 0.09 


1.54 + 0.04 




11.38 


10.0 


0.027 


PKS 1138-26 
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1.07 + 0.04 
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*J2000 co-ordinates from IRAC observations at 3.6/im. 

t850/^m flux densities have been corrected for radio synchrotron e mission contam i nation . 

fThe 850 (im flux density of 12.8 ± 3.3 mjy for PKS 1138-26 from iReuland et al ] l|2004h is the result of a poor quality observat ion. 
References. — (a) DB10 and references therein; (b) Drouart et al. (2012; in prep.) and references therein; (c) 



Archibald et al 



1 20011 ): (d) 



Stevens et al 



(200: 



3). 



Table 2. The first and second order observations with the IRS Long-Low (LL) module for the sample. 



HzRG 


LL First Order 


LL Second Order 




[seconds] 


[seconds] 


3C 239 a 


10 x 2 x 120 


10 x 2 x 120 


4C 23. 56 6 


12 x 2 x 120 




NVSS J180556+633313 a 


20 x 2 x 120 


20 x 2 x 120 


NVSS J171414+501530 a 


20 x 2 x 120 




PMN J0408-2418 c 


15 x 2 x 120 


15 x 2 x 120 


PKS 1138-26 d 


13 x 8 x 120 




4C -00.54 a 


15 x 2 x 120 





Observations. — (a) N. Seymour (PID 50379); (b) G. Fazio (PID 40093); (c) N. Nesvadba (PID 50402); (d) P. Ogle (PID 40460). 
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Table 3. Results from the spectral fitting process with 1 a uncertainties. Columns are as follows: (1) HzRG name; (2) 7.7 ^m PAH luminosity; (3) host galaxy component contribution 
to the best fit in the 24 fim MIPS band; (4) total IR (8-1000 /tin) star-forming luminosity, inferred from 7.7 PAH luminosity; (5) star formation rate, inferred from 7.7 /an PAH 
luminosity; (6) monochromatic rest-frame 6 fim AGN luminosity, (7) mid-IR spectral index of AGN component, (8) 9.7 fim AGN optical depth, (9) specific star formation rate; (10) 
best-fitting iBrandl et al. | j2006h template; (11) x 2 / de 

grees of freedom, v. 
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Mean spectrum 


2.69 ±0.15 


- 


q fil +3-62 

- 1.82 
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shifty 
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u - ui - 0.01 


- 


NGC2623 


162/48 


3C 239 


0.841 ° H 


13 ±3 


1 06 + 110 

- 0.59 


180 1 1 ^ 


4.45 ±0.17 


1 50 + 00 
1 - ou - 0.02 


0.33 ±0.04 


S1 +°- 32 
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Mrk52 
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4.23t° ; 6 6 3 9 
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iooo till 
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u - i0 - 0.15 
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23 if 
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> 1.09 
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36 + 4 
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D - D0 - 3.36 
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i^au _ 580 
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14.40 _ 0AS 


1.41 ±0.07 


0.00 


>0.42 


NGC2623 
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4C -00.54 


2.o 2 ±°; 6 6 I 


21 ±7 
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1 00 + 011 
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0.08l° ; 4 


> 1.19 
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a We take into account the 1 a scatter in the L7.7PAH — ^IR relation. 

Note. — A calibration uncertainty of ~5% is not included for the measurements in columns 2, 4, 5, 6 and 9. 



